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Abstract 



It is shown that, by an appropriate modification of the structure of the interaction po- 
tential, the Breit equation can be incorporated into a set of two compatible manifestly 
covariant wave equations, derived from the general rules of Constraint Theory. The com- 
plementary equation to the covariant Breit type equation determines the evolution law in 
the relative time variable. The interaction potential can be systematically calculated in 
perturbation theory from Feynman diagrams. The normalization condition of the Breit 
wave function is determined. The wave equation is reduced, for general classes of poten- 
tial, to a single Pauli-Schrodinger type equation. As an application of the covariant Breit 
type equation, we exhibit massless pseudoscalar bound state solutions, corresponding to 
a particular class of confining potentials. 

PACS numbers : 03.65.Pm, ll.lO.St, 12.39.Ki. 
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1 Introduction 



Historically, The Breit equation |I| represents the first attempt to describe the relativistic 
dynamics of two interacting fermion systems. It consists in summing the free Dirac 
hamiltonians of the two fermions and adding mutual and, eventually, external potentials. 
This equation, when applied to QED, with one-photon exchange diagram considered in 
the Coulomb gauge, and solved in the linearized approximation, provides the correct 
spectra to order for various bound state problems. 

However, attempts to improve the predictivity of the equation, by developing from it a 
systematic perturbation theory or by solving it exactly, have failed, due to its inability to 
incorporate the whole effects of the interaction hamiltonian of QED 0. Also, this equation 
does not satisfy global charge conjugation symmetry [Q. Finally, the Breit equation, as 
it stands, is not relativistically invariant, although one might consider it valid in the cm. 
frame. 

Nevertheless, despite these drawbacks, the Breit equation has remained popular. The 
main reason for this is due to the fact that it is a differential equation in x-space and thus 
it permits the study of effective local potentials with standard techniques of quantum 
mechanics. Improvements of this equation usually transform it into integral equations in 
momentum space, because of the presence of projection operators 0, ^. 

The purpose of this article is to derive, from relativistic Constraint Theory a 
covariant Breit type equation, where the free part is the sum of individual Dirac hamilto- 
nians. The latter framework ensures the relativistic invariance of the equations describ- 
ing two-particle systems with mutual interactions [^, ^ and establishes the connection 
with Quantum Field Theory and the Bethe-Salpeter equation by means of a Lippmann- 
Schwinger-Quasipotential type equation relating the potential to the off-mass shell scat- 
tering amplitude 0. 

The fact that a covariant Breit type equation can be obtained from Constraint Theory 
was already shown by Crater and Van Alstine [0. In the present paper we show that, 
for general classes of interaction, the covariant Breit type equation is equivalent to the 
Constarint Theory wave equations, provided it is supplemented with a second equation 
which explicitly eliminates the relative energy variable and at the same time ensures 
Poincare invariance of the theory. The potential that appears in the main equation has a 
cm. energy dependence that also ensures the global charge conjugation symmetry of the 
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system. 

The paper is organized as follows. In Sec. 2, we derive the covariant Breit type equa- 
tion from the relativistic Constraint Theory wave equations. In Sec. 3, we determine the 
normalization condition of the wave function. In Sec. 4, we reduce the wave equation, rel- 
ative to a sixteen-component wave function, to a Pauli-Schrodinger type equation, relative 
to a four- component wave function. As an application of the covariant Breit equation, we 
exhibit, in Sec. 5, massless pseudoscalar bound state solutions, corresponding to a par- 
ticular class of confining potentials, involving mainly pseudoscalar and spacelike vector 
interactions. Conclusion follows in Sec. 6. 



2 The covariant Breit equation 

We begin with the Constraint Theory wave equations describing a system of two spin-1/2 
particles composed of a fermion of mass mi and an antifermion of mass m2 , in mutual 
interaction : 

7i.pi - mi^ ^ = - 72.P2 + m2^ , (2.1a) 
-72.^2-^2) § = [ -fi.pi + mi) V ^ . (2.1b) 



Here, is a sixteen-component spinor wave function of rank two and is represented as a 
4x4 matrix : 

^ = ^a,a,{xi,X2) (ai,a2 = 1,...,4) , (2.2) 

where ai{a2) refers to the spinor index of particle 1(2). 71 is the Dirac matrix 7 acting 
in the subspace of the spinor of particle 1 (index ai); it acts on \E' from the left. 72 is the 
Dirac matrix acting in the subspace of the spinor of particle 2 (index 02); it acts on \l/ 
from the right; this is also the case of products of 72 matrices, which act on from the 
right in the reverse order : 



7lM^ = (7^)01/31^/3102 , 72m^ = ^oi/32(7m)/3202 > 
J2^l'l2u^ = ^oi/32(7i^7m)/32Q2 5 = 7^7 7a«,7a/3 



(a = 1,2). (2.3) 



2i 

In Eqs. (2.1) pi and p2 represent the momentum operators of particles 1 and 2, respec- 
tively, is a Poincare invariant potential. 
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The compatibility (integrability) condition of the two equations (2.1) imposes con- 
ditions on the wave function and the potential For the wave function, one finds the 
constraint : 



(Pi - pI) - i^i - ^l) 



^ = , (2.4) 



which allows one to eliminate the relative energy variable in a covariant form. For eigen- 
functions of the total momentum operator P, the solution of Eq. ( |2.4| ) is : 

^ = e-'P-^ ^-t{ml-ml)P.x/{2P^) ^^^t^ ^ (2.5) 
where we have used notations from the following definitions : 

P = Pl+P2, P=^{Pl-P2), M = mi +7712, 

X = '^{xi+X2), X = Xi - X2 . (2.6) 

We also define transverse and longitudinal components of four- vectors with respect to the 
total momentum P : 

qI = Q,-^-^P,. qi = {q.P)P,, P, = P,/Vp^, 



qL = q.P, Pl = VP'. (2.7) 

This decomposition is manifestly covariant. In the cm. frame the transverse components 
reduce to the three spacelike components, while the longitudinal component reduces to 
the timelike component of the corresponding four-vector. (Note that X"^^ = — in the 
cm. frame.) Also notice that, with the definition of the longitudinal components, Pl, 
which is the positive square root of P^, does not change sign for negative energy states 
(under the change P —P)', in this case, it is the longitudinal components, qi, of those 
four-vectors which are independent of P that change sign, since these are linear functions 
of P. 

For the potential, one finds the constraint : 

pI~pI,V^^ = 0, (2.8) 
which means that V is independent of the relative longitudinal coordinate xl '■ 

V = V{x^, Pl, 71, 72) . (2.9) 
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Equations ( |2.5| ) and ( |2.9D show that the internal dynamics of the system is three- 
dimensional, besides the spin degrees of freedom, described by the three-dimensional 
transverse coordinate x^. 

The relationship between the potential V and FejTiman diagrams is summarized by the 



following Lippmann-Schwinger-Quasipotential type |10|, 0, |12[ |13|, |J, |15|, [16|, [l^, i|, |18| 
equation : 



V - T - VGoT 



2Pr. 







T{P,p,p') 



(2.10) 



C(p),C{p') 



where : 

i) T is the off-mass shell fermion-antifermion scattering amplitude; 

ii) C is the constraint 



C{p) = {pI - pI) - {ml - ml) = 2PlPl - {ml - ml 







(2.11) 



in Eq. ( p.lO| ) the external momenta of the amplitude T are submitted to the constraint 
iii) Go is defined as : 

Go(Pi,P2) = ^i(pi) M-V2) Ho , (2.12) 

where Si and 5*2 are the propagators of the two fermions, respectively, in the presence of 
the constraint ( p.ll| ), and Hq is the Klein-Gordon operator, also in the presence of the 



constraint ( |2.11 



^0 = {pl-mD 



{pI - TTT'l] 



O 



p2 1 

/ 2 2\ 



{ml — miy 
4P2 



+ p'' . (2.13) 



c 



In order to obtain the covariant Breit equation, we define the covariant Dirac "hamil- 
tonians" : 



Hi 



-m2l2L - l2Ll2-P2 



(2.14a) 
(2.14b) 



We then multiply Eq. (|2.1a|) by 712, and Eq. ( p.lb|) by 72L , respectively. After subtracting 
the two equations from each other, we obtain the equation : 



(Pl-(Hi+H2))^ 



{PL^{ni^n2)){llLl2LV)^ 



(2.15) 



which can be rewritten as 



Pl{1+1il12lV) - {ni+n2){l-^lLl2LV) 



^ = 



(2.16) 



Addition of the two equations to each other leads to the equation : 



(2.17) 



which can be rewritten as : 



2pl(1-7il72L^) - (Hi-H2)(1 + 7il72L^) 



^ = 



(2.18) 



Upon multiplying this equation by {Tii + 0,2) and noticing that — 7-^2 = "^i — 
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it becomes, after using Eq. ( p.l6|) : 



2pLPL-{ml-ml) (1 + 7il72lV^)^ = 0. 



We now define the Breit wave function \1'b by : 



(2.19) 



= (1 - 7il72l1^)^ • 



(2.20) 



Then, Eq. ( p.l6| ) takes the form : 



Pl{1 + 7il72L^)(1 - 7il72L^)-' - (^i + ^2) 



^1 



(2.21) 



while Eq. (|2.19| ) yields, after factorizing the term (1 + 7il72L^)(1 — 7il72L^) ^ : 



2pl-Pl - {ml - m 



B 



. 



(2.22) 



Equations (|]2TD and (|]2|) are the two wave equations satisfied by the Breit wave 
function b ■ As far as the wave function transformation ( ^.2(J| ) is nonsingular, they are 
equivalent to the initial two wave equations (2.1) of Constraint Theory. 

Equation (|2.21| ) is the obvious generalization of the Breit equation. Its interaction 
dependent part has an explicit cm. energy {Pl) dependence which restores the global 
charge conjugation symmetry that was lacking in the Breit equation. For each solution 
of Eq. ( ^.21| ) with total momentum P , there will correspond, for charge conjugation 
invariant interactions, a charge conjugated solution with momentum —P . 



Equation (|2.22|) determines the relative time evolution law of the wave function, as in 
Eq. (2^), and ensures the relativistic invariance of the theory. While Eq. ( p. 21 ) might 



be considered alone in the cm. frame, Eq. ( ^.221 ) indicates the way of passing to other 
reference frames. 



In the cm. frame, with the standard definitions /3 = 70 and 01=707, Eq. ( p.21|) 
becomes : 

" Po(l + Pip2V){l - f3i(32Vy^ - + ai-p - m2p2 - ct^-p) = ■ (2.23) 

In perturbation theory, V has, in lowest order, according to Eq. ( |2.1CI| ), the structure 

V = --^f/(x^,7i,72) , (2.24) 
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where U is the three- dimensionally reduced form of the propagator of the exchanged 
particle, including the couplings at the vertices. To this order, Eq. ( p. 231) takes the form: 



[ Po - PiP2e{Po)U - (miA + ai.p - ma/^s - c^a-p) ] = . (2.25) 

We notice here, in distinction from the Breit equation, the presence of the energy sign fac- 
tor in front of the potential f/ ; it is this factor which ensures the global charge conjugation 
symmetry of the equation. 

Finally, in the limit when m2 tends to infinity, Eq. (|2.25|) yields the Dirac equation of 
particle 1, with the potential PiU {P2 is replaced by —1 for the antifermion and e(Po) by 
-fl in this limit). 

Equations (2.1), or equivalently ( |2.21| ) and (|2.22| ), were analyzed, in Ref. in the 



nonrelativistic limit, to order 1/c^ , in particular for the electromagnetic interaction case. 
For an arbitrary covariant gauge of the photon propagator, the corresponding hamiltonian 
receives contributions (among others) from quadratic terms generated by the one-photon 
exchange diagram as well as from the two-photon exchange diagrams. However, it turns 
out that in the Coulomb gauge (and also in the Landau gauge to that order) the two- 
photon contribution cancels the quadratic terms arising from the one-photon exchange 
diagram and one then is left with the Breit hamiltonian P, [l^ . This explains why the Breit 
equation in its linearized approximation provides a correct result to order . However, 
in other gauges than the Coulomb and Landau gauges, it is necessary to take into account 
the quadratic terms as well as the two-photon exchange contribution to obtain a correct 
result. 
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3 Normalization condition 



The normalization condition of the wave function \1/ can be determined either from the 
construction of tensor currents of rank two, satisfying two independent conservation laws, 
with respect to Xi and X2 or from the integral equation of the corresponding Green's 
function [|T4|, |T^. One finds for the norm of ip [Eq. (|2.5|) ] the formula (in the cm. frame 



and for local potentials in x^) : 

r ~ r ~ ~ dV 

J £ic Trl ^1 _ v^V + 4710720^0 ■ 



i;} = 2Po, (3.1) 



gp2 

where V satisfies the hermiticity condition : 

= 7io72oV'7io720 • (3.2) 

For energy independent potentials (in the cm. frame) the norm of ip is not positive 
definite for arbitrary V. In order to ensure positivity, it is sufficient that the potential V 
satisfy the inequahty 

^Tr{V^V) < 1 . (3.3) 
In this case one is allowed to make the wave function transformation 



^ (3.4) 



and to reach a representation where the norm for cm. energy independent potentials is 
the free norm. 



In this respect, the parametrization suggested by Crater and Van Alstine |2^, for 
potentials that commute with 7il72l (and hence V'^ = V), 

V = tanhV , (3.5) 

satisfies condition ( p.3| ) and allows one to bring the equations satisfied by ^ [Eq. (|3.4|)] 
into forms analogous to the Dirac equation, where each particle appears as placed in 
the external potential created by the other particle, the latter potential having the same 
tensor nature as potential V of Eq. (|3.5|). 

We shall henceforth adopt the above parametrization (|3.5| ). For more general poten- 
tials that do not commute with 71L72L , the natural extension of parametrization ( p.5|) 
is: 

7il12lV = tanh(7iL72LV^) • (3.6) 



According to Eqs. ( |3.2| ) and ( |3.4| ), we shall introduce the wave function transformation : 

§ = cosh(7ii72Ll^) ^ . (3.7) 



The norm of the new wave function \1/ then becomes (in the cm. frame) : 



1 -U2P2 /^p-7io72oV^_^p7io720^_p7io72oV^_^p-7io72oV^^ 



2Po 

(3.8; 



(The relationship between \1/ and ^p is the same as in Eq. (2.5).) 
Equations (2.1) then take the form : 

( 7i.pi-mi) cosh(7iL72L\^) ^ = (-72-^2 + "^2)7il72L sinh(7ii72L^) ^ , (3.9a) 
{-l2-P2-m2) cosh(7ii72Ll^) ^ = ( 7i-Pi + "^i)7il72L sinh(7iL72Ll^) * . (3.9b) 

In order to determine the normalization condition of the Breit wave function \1'b , we 
first define from V a potential Vb as : 



V, 



B 



1il12lV . 



(3.10) 



(Notice that because of Eq. (|3.2| ) Vg = Vb in the cm. frame.) With this potential, the 
relationship (|2.20|) takes the form : 



B 



cosh Vb 

while the relationship between ^ b and ^ [Eq. ( |3.7| )] is : 



B 



-VB^^ 



The Breit type equation ( p.21|) becomes : 



B 



. 



(3.11) 



(3.12) 



(3.13) 



The normalization condition of ipB (defined from as in Eq. ( |2.5|) ) is, in the cm. 
frame : 



/ 



4 



^2Vb + 2 R 



.2 d 2V; 

0Qp2 



B 



B 



2Po . 



(3.14) 



We therefore end up with three different representations for the two-particle wave 
function. The first one, \l/ [Eqs. (2.1)], corresponds to the framework where Constraint 
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Theory conditions, as well as connection with Quantum Field Theory and the Bethe- 
Salpeter equation, are most easily established. The second one, \I' [Eqs. ( |3.7] ) and (3.9)], 
corresponds to the "canonical" representation, for which the norm, for cm. energy inde- 
pendent potentials, is the free one. This representation is also the one used by Crater and 
Van Alstine |2|]. The third one, [Eqs. ( CT ), ( gl2D , (|3l3D and (|]22D], corresponds 
to the Breit representation. 



4 Resolution of the Breit type equation 

By decomposing the wave function \E' along 2x2 matrix components, Eqs. (3.9) can be 
solved with respect to one of these components and transformed, for the case of potentials 
commuting with ■yiLl2L, into a second order differential equation of the Pauli-Schrodinger 
type |T^. A similar reduction is also possible starting from the Breit type equation ( p.l3| ). 



The relative time dependence of the wave function being determined by Eq. ( p.22| ), 
with a solution of the form ( p. 51) , one decomposes the internal 4x4 matrix wave function 
ipB on the basis of the matrices 1, 7^, 75 and 71,75 by defining 2x2 matrix components : 

4 

i^B = ^/'Sl + 7lV^S2 + 75^53 + 7l75V'S4 = ^ii^Bi ■ (4.1) 

1=1 

We consider the case of potentials that are local in x"^ (but having eventually a cm. 
energy dependence) and that are functions of products of 71 and 72 matrices in equal 
number (general vertex corrections do not satisfy the latter property); then V commutes 
with 71L72L : 

7il72lV = 1^7il72l • (4.2) 
We introduce projection matrices for the above 2x2 component subspaces : 

1 1 

Vi = -(1 + 7il72l) (1 + 715725) , ^2 = -(1 + 7il72l) (1 -715725) , 

Vs = 7(1 - 7il72l) (1 +715725) , 'Pi = 7(1 - 7il72l) (1 - 715725) • 



(4.3) 



They satisfy the relations : 

ViV, = 6i,Vj , V^Tj = 6ijT, (z,j = 1,...,4) . (4.4) 
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(The r's are defined in Eq. (^?T|).) 

Then, the most general (parity and time reversal invariant) potential V [Eqs. ( p.5| ) 
and (|3.10| )] we may consider has the decomposition on the basis ( [4 .31) : 



V = ■ (4-5) 

i=l 

The potentials themselves may still have spin dependences. The spin operators, which 
act in the 2x2 component subspaces, are defined by means of the Pauli-Lubanski oper- 
ators: 

h h 

WiSa = -^ealSfMuP'^Cr'^'' , W2Sa = - ^€aP^^uP'^ (^2'' (e0123 = +1) , 
W^S = W^S = 'l^'P' , Ws = 1^15 + ^^25 . (4.6) 

They also satisfy the relations : 

IllWiSa = ^7?L7l5 , 'y2LW2Sa = ^7^725 • (4.7) 

We introduce the operators : 

f 2 n2 tt/ u/ ( 2 2 Wis-x'^W2s-x'^ 3^ 
w = (^) W,s.W2s , w,2 = (^) ; (4.8) 

then, the potentials a, [Eq. ([4.5D] can be decomposed as : 

ai = A + wB, + wuQ (i = l,...,4), (4.9) 

where the potentials Ai, Bi, Ci are functions of x^"^ and eventually of . 
The projectors (^I3| ) -(f4.5| ) satisfy the simple property : 

exp(^a,P,) = . (4.10) 

^ i=i ^ i=i 

The Breit potential Vb [Eq. ( p. 10 )] has also a decomposition like ( [4.5|) : 



Vb = Y.^B^V^ , (4.11) 

i=l 

with the following relations with the a^'s : 

ttBl = Ctl , (132 = 0.2 , 0_B3 = "0-3 ; 0-54 = —O4, ■ (4-12) 
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The relationship ( ^.121 ) between ips and ip can be rewritten for their 2x2 components 
as well : 

i^Bi = e-«^* ^^Ji (z = l,...,4) . (4.13) 

[ipi is defined from a decomposition of ip as in Eq. ( [4.1| ).] 

The Breit type equation ( |3.13| ) is now easily decomposed into four equations for the 



four components ipBi = 1, • • • ,4) : 

Pl e^^l ^Bl - (mi - 7712) i^B2 + 7-^{WlS - W2s).p ^B3 = 0, (4.14a) 

PLe2«2^^2 - (mi - 7772) V'Bi - 7-^r^s-P^B4 = 0, (4.14b) 
PLe-2«3^^3 _ ^^-54 + T— (W^i5-H^25).pV^m = 0, (4.14c) 
e-2«4 _ m ^bs - nr^^-P = . (4.14d) 

These equations allow one to eliminate the components ipBi, '^B2 and ipBA, in terms of 
ipBs , which is a surviving component in the nonrelativistic limit. Upon defining 

2/7 _ 1 {mj-mlf 2(a,+a,) . .^ 

one finds for ipBi and 7/'b2 the relations : 

Pl^Pbi = e-2(«i + «2 + /^) |_e2«2 (W,s-W,s).p 

{ml - ml) 2a 



[mi — mi) —On^) 
+ -^J^Ws.pe 2«3j^^3, (4.16) 



(7/2^ — 7772) 



p2 

One then obtains two independent equations for ipBs and ipB'i '■ 
MP^e-^^^^l^B, =M'i^B, + (^)V5.pe-2(«i + «2 + /^) 



(1^15-1^25).^ V-Bs • (4.17) 



■2a, 



3 



(1^15 - W,s)-P \ ^B3 , (4.18a) 
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X |-e2«2 {W^s-W2s).p 



+ ^s.p e-2«3| V^Bs . (4.18b) 



Elimination of i^ba leads to the eigenvalue equation for ipB?, '■ 



p2g-2(a3 + 04) _ 



hPr 



2a, 



Ws-p e~2(ai + a2 + h) | ^2ai ^^^ ^ ^ 
{ml — rn^) 



p2 

2 \\-2a 



{W,s-W2s).p t 



hPr 



{W,s - W2s).p e-2(«i + «2 + /^) I e2«2(H^^^ _ W2s).p 

Ws.p e-2«3 } = . (4.19) 



Equation ( |4.19| ) is a second order differential equation for the component i/jbs ■ Usu- 



ally, by wave function transformations one can simplify the structure of the differential 
operators in it. For the general potential ( |4.5| ), the second order differential operator will 
still exhibit a spin dependence. However, for simpler types of potential, the spin depen- 
dence of the second order differential operator also disappears. This is the case of the 
potential composed of general combinations of scalar, pseudoscalar and vector potentials. 
It has the following structure : 

V = V, + 715725^^3 + 71^2 ( gtuV^ + gl^U, + ^T4 ) , (4.20) 
and the decomposition of the potentials [Eq. ([1.5D] along these potentials is given by 



the relations : 

ai = Vi + V2 + V3 + wUi + wi2Ti , 

a2 = Vi + V2 - - wUi - WuTi , 

a3 = Vi - V2 + V3 - wUi - W12T4, , 

= Vi - V2 - V3 + wUi + W12TA . (4.21) 

With potentials of the type ( [4.20 ) and after using the wave function transformation 



= ^Vi-V2 + V3 + 2Ui + 2T^ + h 



14 



X 



{-(2 + 



2Vo - Ua 



T4 



sinh(2?74) 



+ -(1 - wi2)e 



(4.22) 



(wi2 defined in Eq. (|4.8|) ), Eq. ( [4. 191) reduces to a Pauli-Schrodinger type equation, where 
the radial differential operators are those of the Laplace operator. This equation, which 
is also obtained from a wave function transformation in the equation satisfied by ip^ [Eq. 
( [4.13| )], was presented in Ref. [|1^ . 

Equations (4.14) could also have been solved with respect to ipBA instead of ipBs ■ 



5 Zero mass solutions 

As a straightforward application of the covariant Breit equation, with the class of po- 
tentials considered in Sec. 4, we shall exhibit, in this section, a class of solutions which 
correspond to massless pseudoscalar bound states in the limit when the masses of the 
constituent particles tend to zero. 

91/ 

The key observation is that, because of the presence of the kernel e ^ in the nor- 
malization condition (3.14), one is allowed to search for solutions in which some of the 



components ipBi are constants, provided the kernel e is rapidly decreasing at infinity. 

The quantum numbers of the state are determined by those of the components ipBs 
and ipBA , which are the surviving components in the nonrelativistic limit. For the ground 
state they have the quantum numbers s = (for the total spin operator defined in Eq. 
([4.6| )), £ = (for the orbital angular momentum operator) and j = (for the total angular 
momentum operator); these quantum numbers are those of a pseudoscalar state. We shall 
restrict the search by demanding that the components iJjbi and ijjB2 be zero for the ground 
state solution. 

Inspection of Eqs. ( [4.14aj) and ( |4. 14b|) shows that ipBs must be a constant : 



ipB3 = 00 = const. . (5.1) 

(The vanishing of the components iJjbi and 'ijjB2 can then also be checked directly in Eqs. 
(|T|) and (CT)-) 

One is left with the two equations (|4.14c| ) and ( |4.14d| ), which become simple algebraic 



equations : 

PLe~2«30p _ Mt/JBA = , (5.2a) 
15 



PLe 



M0O = . 



(5.2b) 



These equations have a nontrivial solution only if 03 + 04 is a constant : 



03 + 04 = C = const. . 



Then : 



Pl = Me 



M ^° 



c 



(5.3) 

(5.4) 
(5.5) 



We now check the normalizability of the solution thus found. For simplicity, we shall 
consider potentials that are independent of in the cm. frame; the corresponding 
conclusions are not much affected by an eventual smooth P^ dependence of the potentials. 
The normalization condition ( p.l4|) becomes : 



2Pn 



-2a3 



+ e 



2a4 -2C 



2Pn 



8|0opy rf^x e^^as = 2P0 



(5.6) 



where Eqs. ( p.l| ) and (p.3|)-(p3D were used; furthermore, the spin projection must be 
taken in the potentials. Therefore, e~^'^3 must be a rapidly decreasing function when 
|x| ^ 00 , or, equivalently, 03 must be an increasing function of |x| at infinity, indicating 
the confining nature of the potential. 

To have a more explicit representation of the potentials satisfying the above conditions, 
let us consider again the class of potentials composed of general combinations of scalar, 
pseudoscalar and vector potentials [Eqs. ([4.20|) and ( [4.21| )j. Condition (|5.3|) means that 



2{Vi -V2) = C . 



(5.7) 



Thus, the scalar and timelike vector potentials cannot be chosen independently from each 
other. 

The normalizability condition ( ^.6|) implies : 



lim (V3 — 3?74 — T4) = 00 

|x|— >oo 



(5. 
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(The spin projection of the potentials has been taken.) This combination of the pseu- 
doscalar and spacehke vector potentials must therefore be of the confining type. 
The component ipBA [Eqs. (|5.4|) and ( |5.5| )] of the wave function then becomes : 



Equation ( p.5| ) shows that when the masses of the constituent particles vanish, then 
the mass of the bound state also vanishes. [The fact that the right-hand side of Eq. 
(|5.6| ) vanishes in this limit should not lead one to the immediate conclusion that this 
state disappears from the spectrum. It is its coupling to the axial vector current which 
is important on physical grounds, and this coupling involves the relationship of the wave 
function to the Bethe-Salpeter wave function through nonlocal operators, where also Pl 
is involved 0.] 

As far as the potentials do not have singularities at finite distances, the function iIjb4, 
[Eqs. ( ^.4|) and (|5.9| )] does not vanish at finite distances and the corresponding wave 
function ipB does not have nodes; it is then a candidate for the ground state of the 
spectrum. To conclude that this is actually the case necessitates a detailed study of the 
various potentials in all sectors of quantum numbers. Conditions (|5.8|) and (|5.7|) are not 
sufficient to guarantee confinement in general. There are cases of potentials satisfying 
these conditions, for which confinement does not occur in a particular sector of quantum 
numbers or for which some solutions become unnormalizable. However, there are also 
cases for which the above solution is the ground state of the spectrum; in particular, 
when the confining potential is represented by the pseudoscalar potential, Eq. ( [4.19|) 
can easily be analyzed; in this case all solutions other than the one found above remain 
massive in the limit of vanishing constituent masses. 

Equations (4.14) also have solutions for which tpBs = ipBA, = ^ and ipBi and tpB2 are 
nonzero. In this case one finds the solution Pl = \mi — •m2\e~^^'^ ~^ '^2) ^ with (ai + 02) 
equal to a constant. This solution is, however, unphysical, since it belongs to one of the 
unphysical subspaces, where one of the longitudinal momenta, piL or , calculated from 
Eqs. ( |2.6| ) and ( |2.22| ), may become negative 

Finally, the solution found above can also be expressed in the "canonical" representa- 
tion. Taking into account the relationship ( p.l2|) , one finds : 



^54 



e-2(^3-3t/4-T4)^, 



(5.9) 



= (1+ 7^)756 "300 . 



(5.10) 
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The massless pseudoscalar bound state solution found in this section does not of 
course exhaust all possibly existing solutions. Furthermore, several types of mechanism 
may lead to the occurrence of massless pseudoscalar bound states, in connection with the 
spontaneous breakdown of chiral symmetry. Two such mechanisms are : i) the dynamical 
fermion mass generation, due to radiative corrections in the fermion self-energy part |^ ; 
ii) the fall to the center phenomenon, due to short distance singularities . Our solution 
differs from the above two in that it is a direct consequence of the particular confining 
nature of the interaction and therefore hinges on long distance forces, rather than on the 
short distance ones or on the radiative corrections. The solution corresponding to the 
pure pseudoscalar interaction case was studied in detail in Ref. [p3| . 



6 Conclusion 

We have shown that, by an appropriate modification of the structure of the interaction 
potential, the Breit equation can be incorporated into a set of two compatible manifestly 
covariant wave equations, derived from the general rules of Constraint Theory. The 
complementary equation to the covariant Breit type equation determines the evolution 
law of the system in the relative time variable and also determines its relative energy with 
respect to the other variables. Furthermore, in this covariant version of the Breit equation, 
the interaction potential can be systematically calculated in perturbation theory from 
Feynman diagrams by means of a Lippmann-Schwinger-Quasipotential type equation, 
relating it to the off-mass shell scattering amplitude. 

The normalization condition of the Breit wave function indicates the presence of an 
interaction dependent kernel in it, which should be taken into account for consistent 
evaluations of physical quantities, like coupling constants, or for the selection of accept- 
able (normalizable) solutions to the wave equations. In this respect, we exhibited, as a 
straightforward application of the covariant Breit equation, massless pseudoscalar bound 
state solutions, corresponding to a class of confining potentials, essentially composed of 
pseudoscalar and spacelike vector potentials with eventually a particular combination of 
scalar and timelike vector potentials. 

The covariant two-body Breit equation suggests several possibilities for its generaliza- 
tion to the iV-body case {N > 2) or for the incorporation of external potentials. However, 
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one meets here the known difficulty of the "continuum dissolution" problem [El, |2B], which 
prevents the existence of normalizable states. Usually, this difficulty is circumvented by 
the introduction of projection operators, either in the potential ||2^, or in the kinetic 
terms @]. It is not yet known whether some local generalization of the Breit equation 
may avoid the above difficulty. 
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